Topological Defects in Twisted Bundles of Two-Dimensionally Ordered Filaments 
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Twisted assemblies of filaments in ropes, cables and bundles are essential structural elements 
in wide use in macroscopic materials as well as within the cells and tissues of living organisms. 
We develop the unique, non-linear elastic properties of twisted filament bundles that derive from 
generic properties of two-dimensional line-ordered materials. Continuum elasticity reveals a formal 
equivalence between the elastic stresses induced by bundle twist and those induced by the positive 
curvature in thin, elastic sheets. These geometrically-induced stresses can be screened by 5-fold 
disclination defects in lattice packing, and we predict a discrete spectrum elastic energy groundstates 
associated with integer numbers of disclinations in cylindrical bundles. Finally, we show that elastic- 
energy groundstates are extremely sensitive to defect position in the cross-section, with off-center 
disclinations driving the entire bundle to buckle, adopting globally writhing configurations. 
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From the textile fibers in our clothing to the steel 
cables that suspend the Brooklyn bridge, twisted fila- 
ment assemblies are familiar structural elements of ev- 
eryday life. This deceptively simple geometric motif, a 
dense packing of helical filaments with a constant pitch 
and varying radius, underlies several unique mechanical 
properties of ropes and fibers. The twist of a rope or 
yarn transmits tension to lateral compression, binding 
together multiple short strands into single cohesive ele- 
ment of arbitrary length [l|. For the case of wire ropes, 
composed of continuous filaments, helical pre- twist allows 
one to build a material with a high tensile strength while 
maintaining remarkable flexibility from an assembly of 
relatively inextensible elements [2|. It is not surprising 
that Nature has evolved many ways to take advantage of 
the mechanical fidelity of ropes, at much smaller length 
scales, by way of super-molecular assemblies of protein 
filaments in living organisms. Bundles of cytoskeletal 
filaments, f-actin or microtubules, are implicated in a 
variety of cellular processes in eukaryotic life, from cell 
division to mechanosensing [3J. Fibers of extra-cellular 
proteins, like collagen, are ubiquitous mechanical com- 
ponents of plant and animal tissue [4j. Because biolog- 
ical filaments are helically-structured molecules, inter- 
molecular forces between densely-packed filaments gen- 
erate intrinsic torques in bundle assemblies. Indeed, 
there is clear evidence of spontaneous twisting of a va- 
riety biological filament bundle types [5j-|7|], pointing to 
yet another function of twist in these self-organized bun- 
dles. The frustration between two-dimensional packing 
and inter-filament twist generically gives rise to a mech- 
anism of self-limited radial bundle growth |8l-ll0j. sug- 
gesting that spontaneous bundle twist plays a necessary 
role in shaping the structure and mechanics of biological 
filament assemblies. 

While mechanical considerations of ropes and fibers 
date back to as early as Galileo, very little is known 
about the role that the cross-sectional ordering plays in 
determining the global elastic properties of twisted bun- 



dles. In this Letter, we describe a continuum approach 
to the unique, non-linear elasticity of two-dimensionally 
(2D) ordered filament bundles, revealing a fundamental 
coupling between the presence and location of topological 
defects in the lattice packing and the global twist of the 
bundle. We discover a precise equivalence between the 
in-plane stresses created by twisting a bundle and the 
long-range stresses introduced into a 2D elastic sheet by 
forcing it to adopt a spherically curved geometry. Due 
to the accumulation of elastic stresses, we show that a 
sufficiently large and twisted bundle becomes energeti- 
cally unstable to the incorporation of finite number of 
5-fold disclination defects in its cross-section. We pre- 
dict a spectrum of discrete low elastic energy states of 
twisted bundles, associated with a integer numbers of 
disclinations as a function of increasing bundle twist. Fi- 
nally, we show that the global geometry of the low-energy 
states of bundles is extremely sensitive to the location of 
disclinations in the cross-section. Moving the position of 
disclination away from the center of the bundle leads to 
a complex class of structures that buckle, writhing out 
of the plane, ultimately folding the entire bundle into a 
compact helical configuration. 

Consider a bundle of semi-flexible filaments of infinite 
length that is organized into a regular array perpendic- 
ular to mean filament orientation along the z direction. 
Interactions between filaments favor a constant neighbor 
spacing, do, much smaller than the radius of the bundle, 
R, which is approximately cylindrical in shape. The zero 
energy state of these interactions is a hexagonal pack- 
ing of straight filaments. Deviations from the uniform 
spacing as well as filament bending are described by the 
energy [13], 
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Here, A and fj, are elastic constants that penalize com- 
pressive and shear distortions of the array and Uij is elas- 
tic strain tensor of the 2D filament bundle. This strain 



is constructed from displacements of filament positions, 
u(x), in the xy plane, and the projection of the local 
filament orientation, t(x), into the xy plane, 
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The non-linear contribution to strain from filament tilt 
is a direct consequence of the invariance of the elastic 
energy under rotation about in axis in xy plane |ll| , and 
can be intuitively understood to measure distances be- 
tween neighbor filaments perpendicular to the filament 
tangent, that is, the distance of closest approach. The 
final term in eq. (JXJ) describes the bend elasticity of fila- 
ments, with Kb proportional to the stiffness of filaments. 

To study low elastic energy configurations of twisted 
filament bundles whose cross section contains an arbi- 
trary distribution of disclination defects, we focus on con- 
figurations for which the shear and compressive contribu- 
tions to eq. ([l) are constant along filaments. In theses 
structures, u(x) and tj_(x) differ along z only by uniform 
rotations around an axis in the xy plane, and at some ref- 
erence plane (z = 0), the in-plane tilt field is described 
by tj_(x) = A + fl{z x x) where x = is the center of 
the bundle, A is a constant tilt in the plane, and f2 is 
the rate of torsion. Initially, we restrict our analysis to 
A = 0, in which filaments wind around the geometrical 
center of the bundle. 

For fixed fi, configurations are determined by d^ij = 
subject to stress-free boundary conditions at the surface 
of the bundle. As we have reduced the problem to solv- 
ing for the stress <Jij = SE/Suij in a single 2D elastic 
plane, it is convenient to determine the stress tensor 
in terms an Airy stress function, x- By construction, 
o~ij = tiktjfdkdix, is divergenceless, and the Airy stress 
is determined by from "compatibility relation" which re- 
lates \ to the definition of elastic strain, eq. ((2} , 
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where K = 4/i(A + /x)/(A + 2/x) and Vj_ X tj_ = e 
is the 2D curl. Here, s a = ±2tt/6 denotes the topo- 
logical charge of a disclination at position, x a , in the 
cross-section of the bundle where the local bond angle, 
9(x) = 6ijdiUj/2 rotates around singular points. Discli- 
nations represent lattice sites of 5- (positive) and 7-fold 
(negative) symmetry in the otherwise 6-fold symmetric 
packing. The second line of eq. ^ describes the non- 
linear contribution from in-plane tilts to the strain in the 
bundle, deriving from the fact that in 2D ordered chains, 
twisted configurations of tangents necessarily introduce 
gradients in nearest neighbor orientation, lying in a per- 
pendicular plane 13j. For the case of a twisted bundle 
where Vj_ X t± = ±2fl + 0[(flr) 3 ], the second line re- 




FIG. 1: Solid curves show elastic energy of 2D bundle 
groundstates a function of twist times bundle radius, £IR, 
for Kb = 0. The cartoons depict the energy minimizing con- 
figurations of multiple 5-disclinations. The dashed lines show 
the metastable branches of the disclination-free and single- 
disclination bundle energy. 



duces to —K e ff = —eijekididktjti/2 — — 3f2 2 , a uniform 
and negative source of Airy stress in the bundle. 

An analogous compatibility equation describes the 
non-linear elasticity of 2D crystalline membranes, known 
as the Foppl-van Karman theory of thin plates [12]. In 
that case, the coupling between in-plane strain and out- 
of-plane bending generates a source of the Airy stress, 
—K, where K is the Gaussian curvature of the surface. 
Hence, we find the surprising result that twist of fila- 
ment bundle generates precisely the same source of stress 
as stretching a 2D elastic sheet over a positively curved, 
spherical surface of radius R e ff — (VSfl)^ 1 ■ Seung and 
Nelson showed that the far-field stresses created by Gaus- 
sian curvature in 2D crystalline membranes are precisely 
of the form to be "screened" by the presence of point- like 
disclination defects in the membrane [14]. For the case 
of 2D filament bundles, eq. J3j shows that positively 
charged, 5-fold disclinations are able to partially screen 
the elastic stresses of twisted bundles. 

The elastic energy of twisted bundles with an arbitrary 
distribution of disclinations can be computed by solving 
eq. ([3]) for x subject to vanishing normal stress at r = R 
by multipole expansion [15]. These expansions can be 
summed to compute the effective theory for the elastic 
energy, 
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where V is the volume of the bundle and p a is the radial 
distance the ath disclination. The first term describes 
the self-energy of a disclination in the bundle, which no- 
tably, has the same radial dependence as the interaction 
between the disclination and twisted-induced stress de- 



scribed the second term. The third term describes the 
twist dependence of elastic energy in the defect- free state, 
where 7 = K R 2 /K b , that measures the relative impor- 
tance in-plane elastic energy to bend elasticity. The final 
term gives the effective interaction between disclination 
pairs in cylindrical crystals, 
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where |x' - R| 2 = R 2 + p 2 a pj/R 2 - 2p Q ,p |8 cos6», with 
the angular separation between disclinations a and /3. 
The minimum-energy configurations of bundles with 5- 
fold (s a — +2-71-/6) disclinations are shown in FigurcQ]for 
flR < 1.25 and for negligible bending energy, 7 _1 = 
. For these results we find that the elastic energy of a 
filament bundle is unstable to the incorporation a finite 
number of disclinations when the twist of a bundle ex- 
ceeds \fl c R\ = \/2/9 ~ 0.471, at which point the lowest 
energy is accomplished by placing a single 5-fold defect at 
the center. Upon further twist, this single defect struc- 
ture gives way to lower energy states with 2, 3, 4, 5... 
disclinations. Evidently, along with crystals on curved 
surfaces, twisted bundles belong to the unusual class of 
frustrated materials in which defects are necessary com- 
ponents of groundstate configurations. We note that al- 
though 7-fold disclinations are not necessary to screen the 
twist-induced stresses in bundle, we expect these defects 
to appear in low-energy configurations when R 3> dn. 
Just as 5-fold disclinations on large spherical crystals |16| 
are unstable to breaking up into extended strings of al- 
ternating 5- and 7-fold defects of a net positive charge, 
so-called grain boundary scars, we anticipate similar ex- 
tended defect configurations in the elastic-energy ground 
states of large twisted bundles. 

The elastic energy in eq. (j4|) demonstrates an im- 
portant and previously unknown coupling between the 
in-plane stresses created by lattice defects and bundles 
twist: 5-fold disclinations screen the twisted-induced 
stresses and visa versa. Hence, the elastic energy of 
bundles with an excess of 5-fold disclinations (relative 
to 7-fold) can always be lowered by twisting a finite 
amount. The right- or left-handed spontaneous twist 
of a bundle is the analog of the upward or downward 
buckle of a crystalline membrane induced by the pres- 
ence of 5-fold defects in the hexagonal order [14J. The 
degree of "torsional buckling" is sensitive to the rela- 
tive costs of filament bending and in-plane elastic ener- 
gies. We find that optimal degree of twist, 17 1, for a 
bundle with a single centered, 5-fold disclination has the 
form (Oii?) 2 = (3 + 32/7)- 1 . For small bundles of rigid 
filaments, 7 — > 0, the pitch of the bundle sensitive to 
the bending the elastic moduli and insensitive to size, 
2tt/Oi oc y/Ko/Kf,. While for bundles much larger than 
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FIG. 2: The groundstate geometry of elastic energy of eq. 
0, for R = 10do bundles possessing a single disclination 
at radius p. The filled, blue circles show the helical radius 
of the geometrical center of the bundle and the open, red 
circles show the rate of torsion around the z axis. The data 
for each p correspond to disclinations at 6 azimuthal angles: 
0, ±2vr/6, ±4tt/6 and tt. 



this characteristic size, 7 — > 00, the tilt angle of the outer 
filament is predicted by the continuum theory to achieve 
a universal value, VL\R — > \Jl/Z ~ 0.577. 

We now consider filament bundles for which A/0, 
with successive vertical layers are rotated and shifted by 
a constant amount in the xy plane. From the continuum 
theory we compute the mean torque on the filaments due 
to the presence of twist- and defect-induced stress, 
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where aij is the stress in the A = state. Hence, in a 
bundle with a single off-center, 5-fold defect, the elastic 
energy decreases most rapidly by tilt in a direction ir/2 
away from the radial location of the dislocation. When 
A 7^ 0, the center of rotation in the bundle is not the cen- 
tral filament, and eq. §§§ shows that this rotation center, 
for which tj_(x) = 0, migrates from the bundle center 
towards the off-center disclination by a distance, | A|/il, 
demonstrating that the low-energy state of bundle adopts 
a globally writhing and helical configuration. 

To investigate the buckled, helical groundstates of bun- 
dles with off-center disclinations, we introduce a discrete, 
numerical model of cylindrical bundles of fully flexible 
filments with preferred hexagonal order. Bundles are de- 
scribed arrays of filaments at positions, x Q , and filament 
tangents corresponding to, tj_(x) = A + 0(1 x x). The 
elastic energy in a vertical cross-section of height, Az, is 
constructed following the bead-spring model of ref. [lj] , 



AE = \ E (K x « - x ^l - rf o) 2 A^ 
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FIG. 3: Cross-sectional (constant z) and side views of energy 
minimizing bundles for the disclination positions, p < 5do. 
Filament colors correspond to radial distance from center of 
rotation in the xy plane, from yellow r = 0, to blue r — 20do. 



Here, k parameterizes a harmonic restoring force for dis- 
tortions in the packing (x Q — x^j)j_ = x Q — x^j — t(x Q — 
x,g) • t is the projection of the filament spacing perpen- 
dicular to the local tangent direction, thereby, measur- 
ing distances of closest approach for weakly bent fila- 
ments. The interactions are weighted the lengths of ele- 
ments A£ a p = Az(t ■ z)~ l which vary with in-plane tilt 
in the cross-section. To study the groundstates numeri- 
cally, we generate a 5-fold disclination in an unrelaxed, 
infinite hexagonal array. This disclination is offset to a 
distance, p, from the center of the bundle at a given az- 
imuthal angle. Lattice positions within a radial distance, 
r < R = 10do, are triangulated and used as initial in- 
plane coordinates of filaments in the bundle. The elastic 
energy per total filament length is numerically minimized 
by relaxing filament positions, f2 and A. 

The geometry of energy-minimizing configurations is 
plotted in Figure [2] in terms of the rate of torsion, O, 
and the helical radius of the geometric center of the bun- 
dle, r/i, both of which reveal a striking sensitivity of 
global bundle geometry of the bundle to disclination posi- 
tion. The center of helical rotation migrates rapidly from 
Th = when p = to position nearly at the boundary 
of bundle, rn ~ 9.7do when p = bdo, adopting a globally 
helical conformation. Perhaps more surprising, the de- 
gree of spontaneous rotation, Q, steeply increases as the 
disclination position is offset from the bundle center. In- 
deed, the helical pitch of the bundle decreases more than 
3-fold as the disclination is moved from p = to p = 5do. 
Remarkably, the writhing and contorted bundle geome- 
tries shown in Fig. [3] provide a more uniform filament 



spacing than can be achieved in straight bundles. For 
disclinations at radial positions, p > 6do, we find that 
bundle groundstates are self-intersecting, and therefore, 
cannot be modeled without the inclusion of either bend 
elasticity or excluded-volume constraints. 

In summary, we have demonstrated a coupling between 
non-linear elastic stresses that arise from twisting of a 
fillament bundle and orientational defects in the lattice 
packing. Thus, 5-fold disclinations are necessary compo- 
nents of low-energy configurations of sufficiently twisted 
and large bundles of continuous filaments. We expect this 
finding to be particularly relevant to the self-assembled 
bundles and fibers of biological filaments, whose chiral 
molecular structure leads to a generic preference for inter- 
filament twist in compact assemblies. Indeed, cylindri- 
cal collagen fibrils derived from certain tissue types show 
pronounced twist, along with considerable evidence from 
small- angle scattering of localized disorder in their cross- 
sectional organization [6lll7j. In light the results of this, 
we speculate that the unique properties of these self- 
organized collagen fibrils, the well-defined diameter, the 
twisted structure, and their "almost crystalline" order, 
may all find a common origin in generic properties of the 
non-linear elasticity of filament assemblies. 
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